AFRL-ML- WP-TP-2006-4 1 9 

EFFECT  OF  INCLUSION 
MORPHOLOGY  ON  THE 
COEFFICIENT  OF  THERMAL 
EXPANSION  IN  FILLED  EPOXY 
MATRIX  (PREPRINT) 


Kaushik  Mallick,  John  Cronin,  and  Steven  Arzberger 


APRIL  2006 


Approved  for  public  release;  distribution  is  unlimited. 


STINFO  COPY 

This  work  has  been  submitted  to  AIAA  for  publication  in  the  Proceedings  of  the  48th 
Structures,  Structural  Dynamics,  and  Materials  (SDM)  Conference.  One  or  more  of  the 
authors  is  a  U.S.  Government  employee  working  within  the  scope  of  their  Government  job; 
therefore,  the  U.S.  Government  is  joint  owner  of  the  work.  If  published,  AIAA  may  assert 
copyright.  The  Government  has  the  right  to  copy,  distribute,  and  use  the  work.  All  other 
rights  are  reserved  by  the  copyright  owner. 


MATERIALS  AND  MANUFACTURING  DIRECTORATE 
AIR  FORCE  RESEARCH  LABORATORY 
AIR  FORCE  MATERIEL  COMMAND 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OH  45433-7750 


NOTICE  AND  SIGNATURE  PAGE 


Using  Government  drawings,  specifications,  or  other  data  included  in  this  document  for  any 
purpose  other  than  Government  procurement  does  not  in  any  way  obligate  the  U.S.  Government. 
The  fact  that  the  Government  formulated  or  supplied  the  drawings,  specifications,  or  other  data 
does  not  license  the  holder  or  any  other  person  or  corporation;  or  convey  any  rights  or 
permission  to  manufacture,  use,  or  sell  any  patented  invention  that  may  relate  to  them. 

This  report  was  cleared  for  public  release  by  the  Air  Force  Research  Laboratory  Wright  Site 
(AFRL/WS)  Public  Affairs  Office  and  is  available  to  the  general  public,  including  foreign 
nationals.  Copies  may  be  obtained  from  the  Defense  Technical  Information  Center  (DTIC) 
(http://www.dtic.mil). 


AFRL-ML- WP-TP-2006-4 1 9  HAS  BEEN  REVIEWED  AND  IS  APPROVED  FOR  PUBLICATION  IN 
ACCORDANCE  WITH  ASSIGNED  DISTRIBUTION  STATEMENT. 


*//Signature//  //Signature// 

VERNON  BECHEL,  Program  Manager  TIA  BENSON  TOLLE,  Chief 

Structural  Materials  Branch  Structural  Materials  Branch 

Nonmetallic  Materials  Division  Nonmetallic  Materials  Division 


//Signature// 

SHASHI  K.  SHARMA,  Acting  Deputy  Chief 
Nonmetallic  Materials  Division 
Materials  and  Manufacturing  Directorate 


This  report  is  published  in  the  interest  of  scientific  and  technical  information  exchange,  and  its 
publication  does  not  constitute  the  Government’s  approval  or  disapproval  of  its  ideas  or  findings. 

*Disseminated  copies  will  show  “//Signature//”  stamped  or  typed  above  the  signature  blocks. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  searching  existing  data 
sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  this  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a 
collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YY)  2.  REPORT  TYPE 

April  2006 _ Conference  Paper  Preprint _ 

4.  TITLE  AND  SUBTITLE 

EFFECT  OF  INCFUSION  MORPHOFOGY  ON  THE  COEFFICIENT  OF 
THERMAE  EXPANSION  IN  FILLED  EPOXY  MATRIX  (PREPRINT) 


Conference  Paper  Preprint 


6.  AUTHOR(S) 

Kaushik  Mallick,  John  Cronin,  and  Steven  Arzberger 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Composite  Technology  Development,  Inc. 
Lafayette,  CO,  80026 


3.  DATES  COVERED  (From  -  To) 

5a.  CONTRACT  NUMBER 

HQ0006-04-C-7070 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

_ N/A _ 

5d.  PROJECT  NUMBER 

N/A _ 

5e.  TASK  NUMBER 

N/A _ 

5f.  WORK  UNIT  NUMBER 

N/A 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES)  10.  SPONSORING/MONITORING 

Materials  and  Manufacturing  Directorate  AGENCY  ACRONYM(S) 

Air  Force  Research  Laboratory  AFRL-ML-WP _ 

Air  Force  Materiel  Command  ii.  sponsoring/monitoring 

Wright-Patterson  AFB,  OH  45433-7750  agency  report  number(S) 

_ _ ’ _ _ _  AFRL-ML- WP-TP-2006-4 1 9 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

13.  SUPPLEMENTARY  NOTES 

This  work  has  been  submitted  to  AIAA  for  publication  in  the  Proceedings  of  the  48th  Structures,  Structural  Dynamics, 
and  Materials  (SDM)  Conference.  One  or  more  of  the  authors  is  a  U.S.  Government  employee  working  within  the  scope 
of  their  Government  job;  therefore,  the  U.S.  Government  is  joint  owner  of  the  work.  If  published,  AIAA  may  assert 
copyright.  The  Government  has  the  right  to  copy,  distribute,  and  use  the  work.  All  other  rights  are  reserved  by  the 
copyright  owner. 

This  paper  contains  color.  PAO  Case  Number:  AFRL/WS  06-0825,  30  Mar  2006. 

14.  ABSTRACT 

The  paper  presents  material  development  research  performed  at  Composite  Technology  Development,  Inc.  towards  optimization  of  the 
design  for  linerless  composite  cryogenic  tanks.  In  a  cryogenic  composite  tank,  large  thermal  strains  develop  through  the  thickness  of  the 
tank  laminate  due  to  the  mismatch  in  the  coefficient  of  thermal  expansion  of  adjacent  plies  with  different  fiber  orientations.  These  thermal 
strains  are  primarily  caused  by  the  thermal  contraction  of  the  matrix  material.  Excessive  thermal  strains  can  cause  microcracks  in  the 
matrix  and  inter-ply  delamination,  leading  to  leakage  of  the  fluid  contained  by  the  tank.  Reduction  of  the  thermal  expansion  of  the  matrix 
is  viewed  as  an  essential  design  tool  for  optimizing  these  tanks.  Addition  of  inclusions  that  are  much  stiffer  than  the  matrix  is  an  effective 
means  to  reduce  thermal  expansion  of  the  matrix  for  the  composite,  as  long  as  it  doesn't  compromise  the  toughness.  The  paper  presents  an 
analytical  scheme  to  predict  the  effective  thermal  expansion  properties  of  the  filled  matrix  with  embedded  inclusions  and  investigates  the 
effect  of  inclusion  morphology,  shape  and  aspect  ratio  on  the  same.  The  analytical  predictions  are  compared  with  actual  test  results  of 
thermal  expansion  measured  for  the  matrix.  The  inductions  and  trends  are  being  used  to  select  the  best  material  and  to  improve  their 
processing  to  achieve  the  most  optimized  cryogenic  tank  design. _ 

15.  SUBJECT  TERMS 

cryogenic  composite  tank,  thermal  strains,  matrix,  thermal  expansion  properties 


16.  SECURITY  CLASSIFICATION  OF: _  17.  LIMITATION  18.  NUMBER  19a.  NAME  OF  RESPONSIBLE  PERSON  (Monitor) 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  OF  ABSTRACT:  OF  PAGES  Vernon  Bechel 

Unclassified  Unclassified  Unclassified  SAR  16  19b.  TELEPHONE  NUMBER  (Include  Area  Code) 

_ _ _ _J _ _  N/A _ 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39-18 


1 


Effect  of  Inclusion  Morphology  on  the  Coefficient  of  Thermal 
Expansion  in  Filled  Epoxy  Matrix 


Kaushik  Mallick*,  John  Cronin1  and  Steven  Arzberger* 
Composite  Technology  Development,  Inc.,  Lafayette,  Colorado,  80026 


Abstract 

The  paper  presents  material  development  research  performed  at  Composite  Technology 
Development,  Inc.  towards  optimization  of  the  design  for  linerless  composite  cryogenic 
tanks.  Specifically,  it  presents  both  analytical  and  experimental  work  towards  reducing  the 
thermal  expansion  of  epoxy  matrix  for  composite  materials.  In  a  cryogenic  composite  tank, 
large  thermal  strains  develop  through  the  thickness  of  the  tank  laminate  due  to  the 
mismatch  in  the  coefficient  of  thermal  expansion  of  adjacent  plies  with  different  fiber 
orientations.  These  thermal  strains  are  primarily  caused  by  the  thermal  contraction  of  the 
matrix  material.  Excessive  thermal  strains  can  cause  microcracks  in  the  matrix  and  inter-ply 
delamination,  leading  to  leakage  of  the  fluid  contained  by  the  tank.  Therefore,  reduction  of 
the  thermal  expansion  of  the  matrix  is  viewed  as  an  essential  design  tool  for  optimizing  these 
tanks.  Addition  of  inclusions  that  are  much  stiffer  than  the  matrix  is  an  effective  means  to 
reduce  thermal  expansion  of  the  matrix  for  the  composite,  as  long  as  it  doesn’t  compromise 
the  toughness.  The  paper  presents  an  analytical  scheme  to  predict  the  effective  thermal 
expansion  properties  of  the  filled  matrix  with  embedded  inclusions  and  investigates  the 
effect  of  inclusion  morphology,  shape  and  aspect  ratio  on  the  same.  The  analytical 
predictions  are  compared  with  actual  test  results  of  thermal  expansion  measured  for  the 
matrix.  The  inductions  and  trends  from  analytical  computations  are  being  used  to  select  the 
best  material  and  to  improve  their  processing  to  achieve  the  most  optimized  cryogenic  tank 
design. 


Nomenclature 


a  =  aspect  ratio  of  the  inslusion 

CCm  =  coefficient  of  thermal  expansion  of  the  matrix 


OCc  =  coefficient  of  thermal  expansion  of  the  filled  matrix 


Cm  =  elastic  stiffness  tensor  of  the  matrix 

Cf  =  elastic  stiffness  tensor  of  the  inclusion 
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average  strain  in  the  matrix 
volume  fraction  of  the  inclusion 
identity  tensor 
Eshelby  Tensor 
transformation  tensor 

pressure  differential  through  the  tank  laminate 
empirical  constant  in  permeability  testing 
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I.  Motivation 


T  ^inerless  composite  cryogenic  tanks  are  being  considered  by  the  Department  of  Defense  for  applications  that 

can  benefit  from  the  weight,  cost,  and  schedule  advantages  that  these  structures  offer.1  Because  of  their  lower  launch 
mass,  fewer  part  count,  and  quick  manufacturing  turn-around  time,  linerless  composite  tanks  will  be  enabling  for 
launch  vehicles,  especially  in  their  upper  stages,  for  the  ‘Operationally  Responsive  Space  Launch  Initiative’.  It  has 
been  estimated  that  the  use  of  composite  tanks  will  offer  up  to  40  percent  weight  reduction  compared  to  current 
conventional  metallic  tanks,  allowing  increased  payload  capability  and/or  reduced  launch  mass.  Ultra-lightweight 
linerless  composite  tanks  have  also  been  identified  by  NASA  as  enabling  for  future  reusable  launch  vehicles  and  in¬ 
space  propulsion  systems  (Figure  l).2 

Although  composite  materials  hold  much  promise  for  lightweight  tanks,  their  adoption  has  been  slowed  by 
concerns  over  structural  integrity  and  material  compatibility.3  Conventional  composite  materials  typically  become 
brittle  and  suffer  from  microcracking  when 
exposed  to  cryogenic  temperatures,  leading 
to  tank  leakage  or  structural  failure.  The 
use  of  carbon  fiber  based  composites  in 
cryogenic  applications  has  been  explored 
for  the  last  20  years. 3  However,  it  is  only 
very  recently  that  more  novel  and  advanced 
composite  materials  have  been  developed 
and  characterized  well  enough  to  make 
them  viable  alternative  to  metallic 
materials  in  cryogenic  fuel  tanks.  The 
properties  of  these  new  materials  must  be 
explored  further  so  that  they  can  be  used 
with  confidence  and  to  encourage  more  advances  to  be  made  in  future  material  development  efforts. 

One  of  the  technical  challenges  in  placing  composites  in  a  cryogenic  environment  is  to  choose  and  /  or  design  the 
material  to  resist  the  proliferation  of  microcracks  that  may  lead  to  leakage  through  the  wall  of  the  tank  within  the 
tank’s  operating  regime.  Several  studies  in  the  past  have  focused  on  how  composites  microcrack  and  leak  when  they 
are  thermally  cycled  between  cryogenic  and  elevated  temperatures.4'10  Most  of  the  research  conclude  that 
microcracks  are  caused  by  large  thermal  stresses  that  develop  through  the  thickness  of  the  tank  laminate  due  to  the 
mismatch  in  the  coefficient  of  thermal  expansion  (CTE)  of  adjacent  plies  with  different  fiber  orientations.  Since  the 
CTE  of  the  polymer  matrix  is  much  higher  than  that  of  the  unidirectional  fibers,  the  transverse  microcracks  are 
initiated  in  plies  that  want  to  shrink  across  the  fiber  direction  but  are  restrained  by  their  neighboring  plies. 

The  design  of  a  typical  composite  tank  requires  a  laminate  design  that  is  inherently  anisotropic,  resulting  in  the 
CTE  mismatch  between  the  adjacent  plies.  The  feasible  design  optimization  space  that  aims  to  minimize  the  tank 
weight  and  at  the  same  time  mitigate  microcracking  by  manipulating  the  laminate  architecture  is  very  limited.11 
This  leaves  the  cryogenic  tank  designer  with  two  possible  design  routes  to  mitigate  microcracking  in  the  tank  wall, 
both  being  dependent  on  material  development:  1.  improve  the  microcrack  fracture  toughness  and  2.  reduce  the 
thermal  expansion  coefficient  of  the  matrix  material  to  minimize  the  amount  of  residual  stresses  at  cryogenic 
temperature.  The  two  efforts  usually  have  conflicting  demands  on  material  design  and  formulation.  Tremendous 
progress  has  been  made  to  date  at  CTD  in  developing  novel  composite  materials  with  microcrack  resistance 
significantly  better  than  common  off-the-shelf  materials.1  The  current  paper  focuses  on  the  second  option  i.e. 
reduction  of  the  CTE  of  the  matrix  through  judicious  choice  of  inclusions.  Specifically  the  paper  aims  to  achieve 
three  objectives:  1.  identify  the  morphology  of  the  inclusion  that  would  give  the  best  reduction  in  CTE,  2.  correlate 
with  experimental  measurement  of  CTE  of  inclusion  reinforced  matrix  and  3.  optimize  the  reinforcement  volume 
fraction  to  achieve  the  desired  performance. 

II.  Effective  CTE  of  Epoxy  Matrix  with  Embedded  Inclusions 
A.  Model  for  Calculating  the  Effective  CTE 

The  analytical  framework  used  in  the  current  paper  for  the  prediction  of  CTE  is  based  on  the  mean  field  theory. 
The  mean  field  theory  is  a  continuum  mechanics  based  approach  aimed  at  homogenization  of  multiphase 


Figure  1:  Linerless  composite  tanks  will  enable  several  future 
missions  such  as  launch  vehicles  for  Responsive  Space  (left)  and 
reusable  Crew  Exploration  Vehicles  (right) 
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materials.13,14  It  is  used  for  example  in  deducing  their  overall  (“effective”  or  “apparent”)  behavior  (e.g.  stiffness, 
thermal  expansion  and  strength  properties,  heat  conduction  and  related  transport  properties,  electrical  and  magnetic 
properties,  electromechanical  properties,  ...)  from  the  corresponding  material  behavior  of  the  constituents  (and  of 
the  interfaces  between  them)  and  from  the  geometrical  arrangement  of  the  phases.13,13  The  theory  is  based  on  the 
assumption  of  equivalence  of  the  strain  energy  of  the  actual  solid  with  a  disordered  micro  structure  and  that  of  an 
appropriately  defined  effective  continuum.14 

The  model  for  computing  the  effective  CTE  of 
a  matrix  embedded  with  inclusions  is  shown  in 
Figure  2.  The  inclusions  are  assumed  to  be  of  the 
same  size  and  ellipsoidal  shape,  with  an  aspect  ratio 
of  a=a/R  and  transversely  isotropic  properties.  This 
approach  allows  the  computation  to  be  very  flexible 
and  generalized  to  take  into  account  a  wide  variety 
of  reinforcement  morphology  -  (a)  platelet  or  disc 
when  a<l,  (b)  short  fiber  or  whisker  for  a>l  and  (c) 
spherical  when  oc=l.  Takao  and  Taya15  and  Takao16 
had  used  the  Eshelby’s  equivalent  inclusion  method 
to  predict  CTE  of  the  short  fiber/whisker  reinforced 
matrix  material.  This  analytical  approach  is 
extended  in  the  current  effort  to  include  platelet  like 
reinforcement  as  well  as  spherical  reinforcement.  In 
addition  the  analytical  derivations  of  the  effective 
CTE  presented  in  the  proposed  paper  take  into 
account  the  random  distribution  of  the  inclusions  as 

well  as  the  anisotropy  of  both  mechanical  and  thermal  properties  of  the  inclusions  unlike  the  previous  effort. 

The  effective  CTE  of  the  reinforced  matrix  is  derived  as:15,16 

ac=am-  f  | ( TG3T )d,V*e  +  f\ (TG4Aa)dV 

where, 

A  a  =  af  -  am 

OCc  and  OCm  being  the  CTE  of  the  inclusion  and  the  matrix  respectively,  and 

G,  =  {(C, -Cm\s-I)+C,Y‘(C,  -Cj 
G4={(c/-cJs-/)+c/}-'c/ 

c,,c„  are  the  elastic  stiffness  tensors  for  the  inclusion  and  the  matrix,  I  is  an  identity  tensor  and  S  is  the 

Eshelby  tensor  .  The  different  terms  of  the  Eshelby  tensor  S(i,j)  are  provided  in  Appendix  A. 

In  eqn.  (1),  e  defines  the  average  strain  in  the  matrix  due  to  thermal  shrinkage  and  is  obtained  by  solving  the 
following  simultaneous  equations: 

(/  -  /  J  ( TGXT~ 1  )dV  )•  e  +  /  |  {TG2Aa)dV  =  0  (4) 

where, 

Gx  =(S  — I  )G3 

g2  =  (s  ~i)g4  (5) 

In  both  eqns  (1)  and  (4),  T  represents  the  transformation  matrix  provided  in  Appendix  B  and  the  terms  (x'j  represent 
orientation  average  of  the  tensor  v  to  account  for  the  random  orientation  distribution  of  the  inclusions  in  the  matrix.17 

7i  7i  nil 

— 1|  |xsin(a2 ]daxda2da2  (6) 

-7T0  0 
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AT 


Figure  2:  A  model  for  calculating  the  effective  CTE 
of  an  inclusion  reinforced  matrix 
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B.  Calculation  of  Effective  CTE 

Sample  calculations  are  performed  to  determine  the  effective  CTE  of  the  embedded  matrix  for  different  types  of 
inclusions.  Representative  material  properties,  shapes  and  aspect  ratios  assumed  for  the  inclusions  are  listed  in  Table 
1 .  In  the  results  presented,  the  effective  CTE  of  the  embedded  matrix  is  normalized  with  respect  to  the  virgin  matrix 
CTE.  Also  since  the  mean  field  theory  based  derivation  of  the  effective  CTE  assumes  only  a  semi-dilute 
concentration  of  the  inclusions,  the  volume  fraction  of  the  fillers  is  been  limited  to  30%.  This  limit  was  chosen 
arbitrarily,  but  more  specific  limits  will  be  derived  in  the  near  future  based  on  percolation  thresholds  of  inclusion 
connectivity.18, 19 


Table  1:  Elastic  properties  and  thermal  expansion  coefficients  used  in  the  calculations 


Material 

Inclusion 
shape  & 
aspect  ratio 

Elastic  properties 

Thermal  properties 

El 

Ej 

Vlt 

Glt 

aL 

OCt 

GPa 

GPa 

GPa 

ur* 

mm/mm/K 

1  o6 

mm/mm/K 

CTD  500XA 
epoxy  matrix 

- 

3.4 

3.4 

.36 

1.24 

30.0 

30.0 

Nanofibers 

Rod, 

a=300 

400.0 

40.0 

.2 

16.0 

-1.0 

27.0 

Silica 

Spherical,  a=l 

70.0 

70.0 

.2 

29.16 

5.0 

5.0 

Nanoclay 

Disk, 

a=.001 

176.0 

176.0 

.2 

73.3 

3.0 

3.0 

Zirconium 

Tungstate 

Spherical,  a=l 

200.0 

200.0 

.3 

76.92 

-7.2 

-7.2 

Kevlar 

whiskers 

Rod, 

a=300 

61.0 

4.2 

.35 

2.9 

-2.0 

30.0 

PBO  whiskers 

Rod, 

a=300 

270.0 

30.0 

.3 

13.5 

-6.0 

50 

C.  Effect  of  Inclusion  Properties  on  Effective 

Figure  3  compares  the  effective  CTE  of  the 
matrix  material  reinforced  with  several  different 
types  of  inclusions  -  silica  (spherical  particles), 
nanofiber  (rod  shaped),  nanoclay  (disk- shaped), 
Kevlar  whiskers  (rod  shaped)  and  PBO 
whiskers  (rod  shaped).  All  three  of  the  rod 
shaped  inclusions  were  chosen  because  of  their 
negative  axial  CTE  and  strongly  anisotropic 
properties.  The  effective  CTE  of  the  filled 
matrix  was  found  to  be  quite  different  for  these 
three  fibers,  since  the  elastic  as  well  as  thermal 
properties  for  these  fibers  play  a  role  in  the 
effective  CTE  calculation.  For  spherical  and 
isotropic  silica  inclusions,  the  results  indicate  a 
linear  dependency  of  the  effective  CTE  on  the 
volume  fraction  since  it  coincides  with  the 
prediction  of  a  simple  rule  of  mixture.  Of  all  the 
different  inclusions  the  nanofiber  shows  the 
best  promise  in  CTE  reduction. 


CTE 


volume  fraction  (f) 

Figure  3:  Effective  CTE  of  epoxy  matrix  reinforced  by 
different  kinds  of  inclusions. 
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D.  Effect  of  Inclusion  Orientation  on  Effective  CTE 

Figure  4  shows  the  variation  of  effective  CTE  of  the  matrix  material  reinforced  by  nanofibers  which  are  strongly 
anisotropic  in  their  elastic  as  well  as  thermal  properties  (see  Table  1).  The  calculations  are  performed  for  two 
different  inclusion  arrangements  -  parallel  array  and  random.  For  the  former  case,  the  effective  CTE  of  the  material 
in  the  longitudinal  direction  (same  as  the  fiber  axial  arrangement)  is  significantly  lower  but  the  CTE  in  the 
transverse  direction  is  found  to  be  higher 
than  that  of  the  base  matrix  material  for 
small  volume  fraction  of  the  inclusions. 

Such  behavior  is  caused  by  the  axial 
constraint  enforced  by  the  stiff  fibers.  On 
heating  the  composite  (matrix  + 
inclusion),  axial  expansion  of  the  matrix 
is  strongly  inhibited  and  the  resultant 
axial  compression  of  the  matrix 
generates  a  Poisson’s  expansion  in  the 
transverse  direction,  which  more  than 
compensates  for  the  reduction  in  natural 
transverse  thermal  expansion  induced  by 
the  presence  of  the  fibers.20  Similar 
results  have  been  predicted  in  the  past 
literature.15,16  The  MFT  estimates  a  90% 
reduction  in  effective  CTE  along  the 
preferred  orientation  for  a  parallel  array 
and  a  45%  reduction  for  random 
dispersion  of  the  nanofibers,  both  for  a 
5%  volume  fraction  of  inclusions,  which 
far  exceeds  the  expectations  of  a  rule  of 

mixtures.  These  results  also  indicate  that  the  type  of  orientation  of  the  inclusions  which  can  be  influenced  by  the 
processing  conditions,  is  an  important  factor  to  consider  for  estimating  the  effective  thermal  properties. 


Figure  4:  Effect  of  orientation  on  the  CTE  of  nanofiber 
reinforced  epoxy  matrix 


E.  Effect  of  Inclusion  Aspect  Ratio  on  Effective  CTE 

The  expected  reduction  in  the  effective 
CTE  of  the  matrix  reinforced  by  nanofiber  is 
due  to  the  strong  anisotropy  of  the  fibers,  its 
negative  CTE  along  the  axial  direction  of  the 
fibers  as  well  as  its  high  aspect  ratio 
(300-1000).  There  has  been  concern  that  good 
dispersion  of  these  fibers  in  the  matrix  will 
necessitate  high  shear  mixing  which  can 
possibly  reduce  the  aspect  ratio  by  breaking  the 
individual  fibers  during  mixing.  Figure  5  shows 
the  effective  CTE  of  the  reinforced  matrix  for 
different  aspect  ratios  of  nanofibers.  The 
influence  on  the  effective  CTE  is  much  more 
significant  when  the  aspect  ratio  is  reduced 
from  a=  100  to  t2=10  than  it  is  from  <2=1000  to 
a=  100.  The  trend  suggests  that  it  is  extremely 
important  to  control  or  preserve  the  aspect  ratio 
of  the  fibers  during  processing  or  mixing  if  the 
promised  CTE  reduction  is  to  be  realized. 


Figure  5:  Effect  of  aspect  ratio  on  the  effective  CTE  of 
nanofiber  reinforced  epoxy  matrix 
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F.  Effect  of  Inclusion  Shape  on  Effective  CTE 

There  has  been  tremendous  interest  in  Zirconium  Tungstate  (ZrW208),  a  ceramic  with  a  strongly  negative 
coefficient  of  thermal  expansion  (CTE).  In  contrast  to  most  other  ceramics  exhibiting  negative  CTE,  the  CTE  of 
ZrW208  is  isotropic  and  has  a  large  negative  magnitude  (average  CTE  of  -7.2  x  10‘6  mm/mm/K)  over  a  wide  range 
of  temperature  (OK  to  1050  K).21'23  These 
unusual  properties  suggest  the 
incorporation  of  ZrW208  into  a  polymeric 
matrix,  to  create  a  composite  with  very  low 
CTE.  The  objective  of  the  computation 
performed  here  is  to  illustrate  that  the 
inclusion  shape  is  an  important 
consideration,  more  so  than  its  negative 
CTE  in  order  to  achieve  a  matrix  material 
with  low  CTE.  Figure  6  shows  that 
incorporating  spherical  ZrW208  particles 
reduces  the  effective  CTE  only  moderately, 
following  a  rule  of  mixture  as  expected. 

However,  significant  CTE  reduction  can  be 
achieved  if  the  ZrW208  ceramic  can  be 
processed  in  the  shape  of  either  whiskers  or 
platelets.  The  more  pronounced  the  aspect 
ratio,  the  stronger  the  CTE  reduction.  The 
computations  also  show  that  for  random 
dispersion  of  a  given  type  of  inclusion,  the 
most  reduction  in  the  effective  CTE  is 
obtained  with  the  disk  shaped  inclusions. 


0.05 


0.1  0.15  0.2 

volume  fraction  (f) 


0.25 
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Figure  6:  Effect  of  inclusion  shape  on  the  effective  CTE 


III.  Experimental  Results 


A.  Measurement  of  Thermal  Strain 

Specimens  for  measuring  thermal  strain  were  prepared  from  flat  plates  of  CTD-500XA  epoxy  resin  reinforced 
with  different  weight  fractions  of  nanofibers.  32mm  x  32mm  square  test  specimens  were  cut,  prepared  and  strain- 
gaged.  Thermal  expansion  tests  were  performed  using  a  CTD  designed  test  set-up  that  uses  cryogenic  strain  gages 
and  a  comparative  strain  analysis 
technique  developed  by  National 
Institute  of  Standards  and  Technology 
(NIST).  The  NIST  traceable  standards 
used  in  the  test  for  comparison  were 
copper  and  borosilicate  glass  specimens. 

The  thermal  expansion  was  recorded 
from  liquid  Nitrogen  temperature  (77  K) 
to  room  temperature  (293  K).  Each 
sample  was  tested  5  times  and  the  results 
were  averaged. 

Figure  7  compares  the  experimentally 
measured  thermal  strain  of  the  matrix 
embedded  with  different  weight  fractions 
of  nanofibers  and  the  virgin  matrix 
material.  The  results  of  the  tests  show 
that  nanofibers  have  the  potential  to 
significantly  reduce  the  thermal 
expansion  of  the  base  matrix  material. 

There  is  an  overall  reduction  of  24.6%  in 
the  thermal  expansion  of  the  matrix 
using  12%  by  weight  of  nanofiber.  A  summary  of  the  test  results  that  compare  the  reduction  achieved  in  total 
thermal  expansion  at  77K  for  different  concentrations  of  nanofibers  can  be  seen  in  Figure  8. 


Temperature  (K) 

Figure  7:  Thermal  expansion  of  nanofiber  reinforced  CTD- 
500XA  measured  from  77  K  to  293  K. 
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Figure  8:  Reduction  of  thermal  expansion  with  different  weight  fractions  of  nanofiber  inclusion 


CTE  vs.  vol.  fraction  of  VGCF 


B.  Comparison  with  Analytical  Prediction 

Figure  9  compares  the  experimentally  measured  CTE  of  the  nanofiber  reinforced  matrix  with  the  analytical 
predictions.  The  effective  CTE  of  the  matrix 
at  room  temperature  (29 3 K)  was  computed 
from  the  slope  of  the  thermal  expansion  plot 
(Figure  7).  Two  different  aspect  ratios  of  the 
rod-shaped  nanofibers  were  considered: 
a=  10  and  <2=100.  Assuming  that,  the 
experimental  results  agree  fairly  with  the 
analytical  predictions  for  <2=10  for  volume 
fraction  less  than  2.5%.  It  may  be  conjectured 
that  the  aspect  ratio  of  the  nanofibers  is 
significantly  reduced  during  the  processing 
from  its  initial  aspect  ration  of  aspect  ratio  of 
(X  «  300 .  The  nonlinear  trend  in  the  CTE 
reduction  is  encouraging,  with  a  steeper  drop 
for  increasing  concentration  of  nanofibers. 

The  analytical  solution  is  limited  to  dilute 
concentration  of  inclusions  and  is  valid  until 
the  interaction  effect  between  the  nanofibers 
start  to  dominate.  The  ensuing  phenomenon 
needs  to  be  studied  under  the  light  of 
percolation  theory.17 


Figure  9:  CTE  of  nanofiber  reinforced  matrix  at  room 
temperature  (293K) 


IV.  Summary 

Microcracking  in  composite  laminates  subjected  to  thermal  loads  is  an  obstacle  to  the  development  of  linerless 
composite  cryogenic  tanks.  Microcracks  are  a  result  of  large  thermal  strains  that  develop  through  the  thickness  of 
the  tank  laminate  due  to  the  mismatch  in  the  coefficient  of  thermal  expansion  of  adjacent  plies  with  different  fiber 
orientations.  The  paper  presented  both  analytical  and  experimental  work  towards  reducing  the  CTE  of  epoxy  matrix 
for  composite  materials  by  adding  inclusions  that  are  much  stiffer  than  the  matrix.  An  analytical  model  is  first  used 
to  predict  the  reduction  in  the  CTE  of  the  matrix  embedded  with  dilute  concentrations  of  inclusions.  The  effect  of 
the  shape,  aspect  ratio  and  arrangement  of  the  inclusions  on  the  effective  CTE  of  the  matrix  is  investigated 
systematically.  The  trends  in  the  CTE  reduction  derived  from  analytical  calculations  are  used  to  select  the  best  type 
of  inclusions  for  material  development  efforts.  Experiments  are  performed  to  measure  the  thermal  expansion  in 
matrix  embedded  with  nanofibers  from  room  temperature  to  cryogenic  temperature.  Test  results  demonstrate 
significant  reduction  of  the  thermal  expansion  for  moderate  volume  fractions  of  nanofibers.  Finally,  the  effective 
CTE  of  the  matrix  is  compared  with  analytical  predictions  and  inferences  are  made  on  their  agreement. 
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Appendix  A.  Eshelby  Tensors 


The  Eshelby  tensor  S(i,j)  is  defined  by: 


5(U)  = 


1  4dr  -  2 


2(1-0  a2-l 


-  2v_  + 


4  a  -1 


2  +2vm  \g(a) 


5(2,2)  =  5(3,3)  =  ——— 


4(l-vm)  2{a2- 


4a2  -13 


-2 vm  \g{a) 


5(1,2)  =  5(1,3)  =  — - 


2(1-0  1  -a 


-  +  2v  + 


2a2  +1 


-2vm  g(ar) 


5(2,1)  =  5(3,1)  = 


1  2a2  2a2  +1 


4(1  - v  )  1-dr  a  - 1 


+  2vm  h(ar) 


S(23)  =  S(3,2)=— !- 


4(l-vm)  2«'-l  4(1  -  a 


4 a2  -1 


2  \ +  2vm  |s(«) 


5  (4,4)  =  5  (5,5)  =  ——— 


2d -vm)  1-flr 


■  -  2v.„  +  - 


1  2a2 +4 


■5(6,6)  =  -— - 


4 a2  -7 


2  or  -1 


-2vm  g(ar) 


+  2vm  Uar) 


v  ,  /  2(1 -vm)  \  2\cc2 -1)  i_4(dr2  -i)  -r  'j 

all  other  components  of  S{i ,  j)  being  zero.  The  function  is  defined  as 

g{a)  =  - — ^  [dr(dr  - 1)1/2  -  cosh-1  (or)  for  a  >  1  ( 


s(ff)  = 

g{a)  = 


[a2  - 1)3 

(l-a2)3 


for  OC  >  1  (whiskers/fibers)  (A2) 


(df)  —  d"(l  —  a1 )  for  (X  <  1  (platelets/disks) 


For  d!  =  1  (spherical  inclusions),  the  components  of  S{i,  j )  can  be  calculated  by  limit  investigation  for  i,  j  <  3  as 


5(i,i')  = 


7  -  5v_ 


15(1- vj 


5 O',  j)  -  - 


1  -  5v„ 


for  i  i1  j 


15(1-0 


,V(i. i )  —  7  fori,  7  >3 

15(1-0 
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Appendix  B.  Transformation  Matrix 

The  transformation  matrix  T  relates  the  local  coordinate  system  of  the  inclusion  ( X ^  )  to  the  global 

coordinate  system  (xj ,  X2 ,  X3  )  and  is  given  by:17 


Mi  12 

T  T 

1  21  1  22 


r„  = 


L 12 


T  = 

1  21 


L  22 


where, 
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sin  a,  cos  +  cos  ax  cos  a2  sin  a3 


rl3  =  sin  a2  sin  a3 


r2l  =  -  cos  ax  sin  a3  -  sin  ax  cos  a2  cos  a3 

r21  =  _  Sin  a\  Sin  a3  +  C0S  a\  C0S  a2  C0S  a3 


(B2) 


r23  =  sin  a2  cos  a3 
r3l  =  sin  ax  sin  a2 
V32  -  -COSOTj  sin6f2 


r33  =  cos  a2 
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